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Values  of  Diffusion  Coefficients  Deduced  From 
the  Closing  Times  of  Helicopter-Produced 

Clearings  in  Fog 

1.  INTRODUCTION 

Elliott1  indicated  methods  whereby  values  of  the  mixing  coefficients  in  stable 
atmospheric  situations  might  be  estimated  from  the  closing  times  of  clearings 
created  in  fog  layers  by  artificial  means.  He  considered  the  case  in  which  a slot 
of  clearing,  of  sufficient  length  and  depth  to  be  considered  infinite,  had  been 
created  in  a fog  layer,  and  presented  equations  describing  the  "closing  in"  charac- 
teristics of  the  slot  as  a function  of  the  initial  slot  width  and  time.  He  also  pointed 
out  that  similar  equations  could  be  written  for  a circular  hole  of  clearing  which 
could  be  treated  mathematically  as  an  infinite  cylinder. 

Clearing  data  that  can  be  employed  to  estimate  the  values  of  diffusion  coeffi- 
cients in  the  manner  of  Elliott's  suggestion  have  been  acquired  by  Plank,  Spatola, 

2 3 

and  Hicks.  ’ The  clearings  were  created  by  helicopter  downwash  during  hover 
experiments  performed  5 miles  north  of  Lewisburg,  West  Virginia,  in  September 

(Received  for  publication  11  January  1977) 

1.  Elliott,  W.P. , (1970)  The  role  of  diffusion  in  closing  artificially-produced 

holes  in  clouds,  J.  Appl.  Meteor.  9:900-902. 

2.  Plank,  V.  G. , Spatola,  A.  A.,  and  Hicks,  J.R.  (1970)  Fog  Modification  by  Use 

of  Helicopters.  AFCRL-70-0593;  also  ECOM  R & D Tech.  Rept.  5339, 

58  October  1970,  154  on. 

3.  Plank,  V.  G. , Spatola,  A. A.,  and  Hicks,  J.R.  (1971)  Summary  results  of  the 

Lewisburg  fog-clearing  program,  J.  Appl.  Meteor.  10:763-779. 
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1969.  The  geometric  and  thermodynamic  properties  and  observed  closing  times  of 
certain  of  these  clearings  are  summarized  herein,  and  the  modifications  to  Elliott's 
problem-treatment  method  required  to  estimate  the  diffusion  coefficients  for  the 
particular  Lewisburg  situations  are  discussed.  The  estimation  methods  are  indi- 
cated and  the  values  presented,  with  qualifying  comments  regarding  their 
significance. 

2.  THE  CLEARING  DATA 

Nine  helicopter -hover  experiments  were  conducted  during  the  Lewisburg  pro- 
gram, and  they  yielded  data  that  were  suitable  for  estimating  the  mixing  coefficients. 
The  dates  and  times  of  these  experiments  are  listed  in  Table  1 together  with  infor- 
mation about  the  geometric  sizes  and  thermodynamic  properties  of  the  clearings  at 
the  surface  and  at  fog-top  levels.  The  measurements  were  obtained  under^steady- 
state  conditions  of  clearing,  as  described  by  Plank,  Spatola,  and  Hicks,  ’ or 
they  were  derived  by  methods  indicated  in  the  table  caption.  The  ESTs  listed  are 
for  when  the  helicopters  ceased  hovering  and  when  aircraft  observations  of  the 
closing  of  the  clearings  first  began.  The  closing  times  of  the  clearings  are  shown 
in  the  last  column  of  the  table. 

The  clearings  had  truncated,  conical  shapes  and  the  air  temperature  within  the 
clearings  was  warmer  than  that  of  the  surrounding  fog  due  to  the  adiabatic  warming 
of  the  helicopter  downwash  air  and  the  heat  of  the  engine  exhaust  that  was  added  to 
the  air.  The  water  vapor  content  of  the  clearings,  although  less  than  saturated, 
was  likewise  greater,  in  all  but  one  case,  than  the  water  mass  content  (vapor  plus 
liquid)  of  the  surroundings.  This  was  the  result  of  the  particular  fog  situations  at 
Lewisburg  in  which  the  clear  air  above  the  fog,  which  was  transported  downward 
into  the  fog  by  the  helicopter  rotors  during  the  clearing  process,  had  a larger  water 
mass  content  than  that  of  the  fog  air  itself. 

3.  MODIFICATIONS  TO  THE  PROBLEM  TREATMENT  OF  ELLIOTT 

The  equations  of  Elliott  cannot  be  applied  directly  to  estimate  the  diffusion 
coefficients  for  the  Lewisburg  clearings  because  (a)  they  pertain  to  a slot,  or  swath, 
of  clearing,  whereas  the  Lewisburg  clearings  had  conical  shape,  and  (b)  they 
describe  the  diffusion  of  water  substance  in  isothermal  situations,  whereas  the 
Table  1 data  reveal  that  the  Lewisburg  clearings  were  distinctly  non-isothermal  in 
character.  These  differences  require  treatment  modifications  as  indicated  in  the 
following  sections. 
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Table  1.  Dates,  Times 


t.  Dim  SION  Kyi  VTiONS  OK  ft  YTKH  >1  IM  \M  K \M)  IIKVT  KOK  \\ 
IM'IMTK  KYKINDKIl 


It  is  assumed,  for  first  estimation  purposes,  that  the  conical-shaped  clearings 
of  the  Lewisburg  experiments  can  be  approximated  by  the  theoretical  diffusion  equa- 
tion for  an  erect,  right  cylinder  of  sufficient  length,  that  is,  fog  depth,  to  be  con- 
sidered infinite.  (The  details  of  how  the  cones  were  approximated  by  cylinders  is 
described  in  Section  7. ) The  governing  equation  that  describes  the  diffusion  of 
water  substance  into  or  out  of  the  interior  of  the  cylinder  and  from  or  to  the  bound- 


ary (in  the  temporarily  presumed  absence  of  any  second  interacting  diffusion  proc- 

4 

ess)  is,  in  cylindrical  coordinates. 


9C(r,  t) 

at 


K, 


a2C(r,  t) 

— 2 


9r 


I 3C(r,  t) 
r 9r 


(1) 


where  r is  radial  distance  outward  from  the  centerline  of  the  cylinder,  t is  time 
measured  from  some  initial  time  t = 0,  K^,  is  the  eddy  diffusion  coefficient,  and 


C = pq  + M 


(2) 


is  the  "water  substance  parameter,"  where  p is  the  air  density,  q is  the  specific 


humidity  of  the  water  vapor  in  the  air,  and  M is  the  liquid  water  content  of  the  fog, 

1 


when  present.  Equation  1 is  analogous  to  that  employed  by  Elliott,  and  its  appli- 
cation to  the  fog  problem  is  based  on  the  assumptions  that  the  coordinate  system 
moves  with  the  cloud,  that  horizontal  wind  shear  effects  can  be  ignored, . that  the 
eddy  diffusion  of  water  vapor  and  liquid  water  in  droplet  form  are  both  governed  by 


the  same  constant  coefficient  K^,,  and  that  evaporation  and  condensation  occur 


instantaneously. 

If  differences  of  temperature  (T)  exist  initially  between  the  interior  of  the 
cylinder  and  its  boundary,  and  if  we  assume  temporarily  that  no  second  diffusion 
process  is  operating  that  would  modify  the  heat  flow  by  interaction,  then  the  tem- 
poral and  spatial  changes  of  temperature  within  the  cylinder  will  be  governed  by  the 
equation 


9T(r,  t)  ^ 
5t K 


-[■ 


9 T(r,  t) 


9r 


I 9T(r,  t) 
r 9r 


(3) 


4.  Carslaw,  H.  S. , and  Jaeger,  J.C.  (1959)  Conduction  of  Heat  in  .Solids. 
Clarendon  Press,  Second  Edition. 


Oxford, 
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which  is  of  type  similar  to  Eq.  (1),  where  KH  is  the  eddy  coefficient  of  thermo- 
metric conductivity. 

The  solutions  of  Eqs.  (1)  and  (3)  for  any  given  set  of  initial  conditions,  identi- 
fied by  the  subscript  "0,"  and  boundary  conditions,  identified  by  the  subscript  " 1,  " 
are,  respectively^ 


, * J (rff  ) 2 

1 - - Y — j tz x X exp  (-  K„  a t) 

a Zv  J.laa'  ) r C n 

, n 1 n 

n=  1 


(4) 


which  satisfies  the  initial  conditions  C = Co  at  a > r > 0,  at  t - 0,  and  boundary 
conditions  C = at  r = a,  for  t > 0,  and 


1 - 


n=  1 


a-TTld-)  Xexp(-KHa2t) 
n l n 


(5) 


which  satisfies  analogous  initial  and  boundary  conditions  on  T.  In  these  equations, 

.1  is  the  Bessel  function  of  zero  order  of  the  first  kind,  J . is  the  Bessel  function 
o 1 

of  the  first  order,  which  fulfills  the  boundary  conditions,  and  a^,  0-2,  

are  the  positive  roots  of  (oa)  = 0. 


5.  NOVISOTIIERM \l„  SIM1I.T  WKOliS  1)11  I I SION  WITHOI  T COI  IM  IM. 


Henry5  and  Frank-Kamenetskii6  have  demonstrated  that  if  two  diffusing  sub- 
stances are  involved  in  a given  problem  (having  its  particular  geometric,  boundary, 
and  initial  conditions),  the  problem  may  be  regarded  as  one  in  which  two  time- 
waves  of  diffusion  occur  for  each  of  the  substances.  There  is  a "primary  wave,  " 
which  moves  independent  of  the  existence  of  the  other  substance  and  is  governed  by 
equations  of  the  type  (1)  and  (3)  herein,  and  there  is  a "coupled  wave,  " which  moves 
dependent  on  the  mutual  interaction  influences  of  the  two  substances. 

If  we  apply  the  Henry5  findings  to  the  present  problem,  in  which  we  have  a 
cylinder  of  clearing,  and  if  we  assume  (a)  that  the  boundary  conditions  within  the 

5.  Henry,  B.A.  (1939)  Diffusion  by  absorbing  media,  Proc.  Royal  Soc.  171: 

215-241. 

6.  Frank-Kamenetskii,  D.A.  (1955)  Diffusion  and  Heat  Exchange  in  Chemical 

Kinetics,  Princeton  University  Press,  Princeton,  N.J.  (translated  by 
N.  Thon). 
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environment  are  invariant  with  time,  (b)  that  the  coupling  interaction  between  the 

$ 

diffusion  of  heat  and  the  diffusion  of  water  substance  is  negligible,  and  (c)  that  the 
two  eddy  coefficients  of  diffusion  are  equal,  constant  with  time,  and  independent  of 
the  magnitude  changes  of  either  the  temperature  or  water  vapor  content  of  the  air. 


Kc  = kh  = Ke  , (6) 

then  we  may  conclude  that  the  initial,  boundary,  and  variable  parameters  will,  at 
all  times  during  the  diffusion  processes,  be  related  as 


C - C T - T 
o _ _o 

C - C " T - T. 


This  relationship  implies,  in  essence,  that  the  fractional  time  changes  of  tem- 
perature and  water-substance  within  the  cylinder,  at  any  point,  will  proceed  at  the 
same  rates,  proportional  to  their  initial  differences. 


Figure  1.  Solutions  of  Equation  (4)  or  (5)  for  Two  Ranges  of  the  Ratio  (C  - C) / 
(CQ  - Cj)or  (T0  - T)/(T0  - Tj).  The  isolines  give  the  values  of  the  non-aimen- 
sional  quantity  Ket/a2.  Presented  with  the  permission  of  the  Clarendon  Press, 
Oxford,  England,  publisher  of  the  book  "Conduction  of  Heat  in  Solids,"  by  H.S. 
Carslaw  and  .T.  C.  Jaeger,  1959. 

The  coupling  between  the  water  substance  and  heat  in  the  fog  problem  involves  the 
latent  heat  released  by  the  conversion  of  water  vapor  into  liquid  water  during  the 
mixing  processes.  The  amount  of  this  heat  release  is  small  and  can  be  ignored, 
to  a first  approximation.  For  example,  the  condensation  of  fog  liquid  water  in 
amounts  of  the  order  of  0.  10  gm/m*  in  the  temperature  range  5°  to  20°C,  will 
only  cause  temperature  changes  of  about  0.  1 C°. 


Plots  of  Eqs  (4)  and  (5)  that  incorporate  assumptions  (6)  and  (7)  are  shown  in 

2 

Figure  1.  The  isolines  specify  the  values  of  the  non-dimensional  parameter  Ket/a 
which  correspond  to  particular  values  of  the  ratio  quantities  (CQ  - C)/(Cq  - Cj)  or 
(T  - T)/(Tq  - Tj)  and  to  particular  radial  distances  r/a  from  the  centerline  axis 
of  the  cylinder.  The  diagram  is  a conversion  of  Figure  24  of  Carslaw  and  Jaeger 
and  is  analogous  to  the  plot  discussed  by  Elliott.  * 


6.  SATURATION  ASSUMPTIONS,  CONCERNING  THE  THRESHOLD  OF 
CONDENSATION  AT  ANY  GIVEN  RADIUS  AND  TIME  WITHIN  THE 
CYLINDER 

It  is  assumed  that  the  fog— no  fog  boundary  within  the  cylinder  will  move  rad- 
ially inward  with  time,  in  accordance  with  the  observations  of  the  Lewisburg  exper- 
iments, and  that  the  air  at  any  particular  point  on  this  boundary,  at  the  radius  rg 

and  time  t , when  condensation  first  occurs,  will  be  saturated  but  devoid  of  liquid 
s’ 

water,  such  that 

M = 0 (8) 

and,  from  Eq.  (2), 

c = cg  = Psqs  • <9) 

If  the  pressure  in  the  cylinder  at  any  horizontal  level  is  considered  to  be  the 
same  as  that  of  the  environment,  the  air  density  at  saturation,  pg,  is  given  by  the 
equation  of  state  as 


ps  = R(1  + 0.  6078  q )T 
s s 


(10) 


where  R is  the  gas  constant  for  dry  air  and  Tg  is  the  absolute  temperature,  in  °K, 
at  which  saturation  occurs.  The  saturation  specific  humidity  and  the  saturation 

7 

temperature  are  related  by  the  equation  of  Tetens, 


7.  Tetens,  O.  (1930)  Uber  einige  meteo.  ologische  Begrtffe,  Z.  Geophys. 
6:297-309.  
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where  P^  is  specified  in  millibars  and  the  constants  are  those  of  Murray. 


If  the  initial  and  boundary  conditions  are  specified  for  any  given  fog  situation 
for  which  the  cylinder  approximation  is  presumed  to  apply,  Eqs  (7),  (9),  (10),  and 
(11)  may  be  solved  simultaneously  to  determine  the  particular  values  of  Tg,  q 
and  Cg  that  pertain  to  the  "threshold  of  condensation"  condition  within  the  cylinder. 
The  Cg  and  Tg  values  may  then  be  introduced  in  the  ratio  quantities  on  the  left  of 
diffusion  Eqs  (4)  and  (5).  From  these  diffusion  equations,  thus  evaluated,  the 
radius-time  characteristics  of  the  diffusion  within  the  cylinder  can  be  ascertained, 
as  is  illustrated  in  Figure  1. 

7.  DETERMINATION  OE  THE  DIFFUSION  COEFFICIENTS 

It  was  assumed  that  the  conical-shaped  clearings  of  the  Lewisburg  experiments 
could  be  approximated  by  cylinders  having  a radius  a equal  to  the  average  of  the 
surface  and  fog-top  level  radii  that  were  observed.  These  a values  are  listed  in 
Table  2,  together  with  the  mean  values  of  the  initial  and  environmental  parameters 
To,  qQ,  Pj,  Tj,  q^,  and  M^.  These  means,  except  for  M^,  are  the  vertical  aver- 
ages of  the  measured  surface  and  fog-top-level  values  of  Table  1.  Fog  liquid  water 
was  measured  only  at  the  surface  level,  hence  the  Mj  values  are  those  of  the 
surface. 

The  closing  times  of  the  clearings,  t , are  also  listed  in  Table  2:  it  should  be 
noted  that  these  correspond  to  the  second  of  the  observed  times  of  Table  1.  (The 
first  of  the  times  specified  in  Table  1 is  the  approximate  time  required  for  the 
initial  whisps  of  fog  to  diffuse  inward  to  the  center  of  the  clearings  from  the  sur- 
roundings; the  second  is  the  time  required  for  the  vertical  visibility  at  the  center 
of  the  former  cleared  zones  to  become  approximately  the  same  as  that  of  the  sur- 
roundings. It  is  assumed  herein  that  the  latter  times  are  the  most  appropriate  for 
estimating  the  diffusion  coefficients  for  the  layer-mean  conditions  of  clearing  under 
the  cylinder  approximation. ) 

8.  Murray,  F.  W.  (1967)  On  the  computation  of  saturation  vapor  pressure, 

J.  Appl.  Meteor.  6:203-204. 
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With  these  assumptions  about  the  correspondence  of  the  observational  data  and 
theory,  the  Tg,  qg,  and  Cg  values  were  computed  for  each  of  the  nine  experiments, 
following  the  method  indicated  previously.  Absolute  accuracy  of  the  measured 
parameters  was  assumed  in  the  computations,  and  it  might  be  mentioned  that  a high 
degree  of  computational  exactitude  was  required.  The  resultant  values  are  given 
in  Table  2. 

Knowledge  of  these  values  permitted  determination  of  the  ratios  (C  - C )/ 

o s ' 

(C  - C.)  or  (T  - T )/T  - T.),  and  this,  in  turn,  from  Eqs  (4)  and  (5),  as  shown 

OX  U b O 1 n 

plotted  in  Figure  1,  gave  the  (Kgt/a  >r_p  values  that  pertained  to  the  centerline  of 
the  cylinder,  where  r = r/a  = 0.  Finally,  the  diffusion  coefficients, 

Ke  = (Ket/a2)r=0Xa2/tc  , (12) 

2 

were  computed  from  the  thus -established  values  of  (K  t/a  ) and  from  the 

e r = 0 

observed  a and  t„  values.  The  K values  are  listed  in  Table  2. 
c e 

8.  DISCUSSION 

It  is  seen  that  the  K values  for  the  Lewisburg  experiments  ranged  from 
5 5 e 2 

0.  7 X 10  to  1.9  X 10  cm  /sec.  The  values  were  reasonably  consistent  for  exper- 
iments performed  on  the  same  days,  but  there  was  little  correlation  with  the  visual- 
appearance -state  of  the  fog  top,  or  with  any  of  the  other  meteorological  or  geometric 
parameters  identified  in  Tables  1 and  2. 

These  K values  are  larger,  by  about  a factor  of  two,  than  the  profile  values 
c 9 

reported  by  Vorontsov  and  Selitskaia  for  250-meter-thick  fog  in  the  Dickson 
$ 

Island  region  during  landward  wind  conditions,  and  are  likewise  larger  than  the 
values  for  sea  fog  reported  by  Vorontsov.  * The  values  correspond  approximately 
with  the  typical  values  for  stratus  clouds  (0.44  X 105  cm2/sec)  and  altostratus 
clouds  (1.  1 X 10  cm  /sec),  as  reported  by  Smith,  Chien,  and  MacCready.  ^ 


A small  island  off  the  NW  coast  of  Siberia  at  the  mouth  of  the  Yenisei  River. 

9.  Vorontsov,  P.  A.,  and  Selitskaia,  V.  I.  (1959)  The  vertical  structure  of 
summer  fog  in  the  Dickson  Island  region,  Leningrad.  Arkticheskii  i 
Antarkticheskii  Institut,  Trudy  228(1):87-5TT 

10.  Vorontsov,  P.A.  (1959)  Aerological  characteristics  of  sea  fog,  Leningrad. 

Arkticheskii  i Antarkticheskii  Nauchno-Issledovatel'skii  Institut  Trudv 

~Z2B(l):3fl-54; — 1 1 

11.  Smith,  T.B.,  Chien,  C. -W.,  and  MacCready,  P.  B.  (1969)  Study  of  the 

Engineering  of  Cloud  Seeding,  Final  Report,  Contract  N00189-68-C-0100, 
Meteorology  Research,  tnc.7  Altadena,  Calif. 
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There  is  little  question  but  that  the  tabulated  values  herein  are  the  "effective 
coefficients'  that  pertained  to  the  closing-in  of  the  helicopter  clearings.  But  there 
is  a question  as  to  whether  the  values  are  characteristic  of  the  ambient  fog 
surroundings. 

There  are  several  obvious  reasons  why  the  computed  values  might  be  too 
large,  relative  to  the  true  ambient  values.  But  there  is  also  a counter  reason  for 
arguing  that  they  might  be  proper  or  too  small.  The  actual  truth  is  unknown,  of 
course,  and  we  can  only  surmise  the  possibilities. 

The  tabulated  values  could  readily  be  overestimates  of  the  true  ambient  values 
because,  at  the  initial  times  and  during  the  closing-in  process,  residual,  helicop- 
ter-induced turbulence  could  have  been  present  in  the  clearings  that  was  not  taken 
into  account  in  the  theoretical  treatment.  Moreover,  since  the  clearings  were 
initially  warmer  than  the  surroundings,  solenoidal-type  circulations  could  have 
existed  across  the  clearing  boundaries,  which  is  another  form  of  turbulence  not 
considered.  Wind  shear  effects  could  also  have  operated  to  close  the  clearings 
more  rapidly  than  in  the  homogeneous  wind  field  situation  assumed  in  the  theoret- 
ical treatment,  again  leading  to  an  overestimate  of  the  coefficient  values.  On  the 
other  hand,  though,  the  equations  employed  herein,  which  presume  the  cylinder 
boundary  to  be  fixed  at  the  radius  a,  with  invariant  environmental  conditions  pre- 
vailing at  and  beyond  this  boundary,  are  not  really  good  approximations  of  the  true 
situation.  In  the  real  atmosphere,  as  in  the  Lewisburg  experiments,  it  would  seem 
that  the  eddy  diffusion  of  water  substance  and  heat  across  the  clearing  walls,  which 
existed  initially  at  t = 0,  would  proceed  in  both  radial  directions  with  time,  both 
toward  center,  toward  decreasing  r,  and  also  outward,  toward  r > a.  If  this  were 
true,  it  would  imply  that  the  values  of  the  table  might  be  underestimates  of  the 
actual,  by  some  unknown  factor  related  to  this  boundary  specification  problem. 

In  any  event,  it  may  be  stated  that  the  clearings  of  the  Lewisburg  experiments 
did  close-in  very  quickly,  which  is  prima  facie  evidence  of  large  turbulent  diffusion. 
But  whether  the  method  of  Elliott1  provides  the  proper  coefficient  values  for  the 
diffusion  remains  uncertain.  Comparison  experiments  are  needed  in  which  the 
coefficients  deduced  from  hole-closing  are  related  directly  to  those  of  conventional 
measurement. 


SI  \1\I\R\  ll.l.l  STR  \TIOYS 

The  salient  features  of  the  diffusion  situations  discussed  in  this  paper  are 
illustrated  in  figures  2 and  3.  The  upper  diagram  of  Figure  3 shows  how  the 
closure  times  of  a cylindrical  clearing  having  the  mean  properties  stated  in  the 
caption  will  vary  as  a function  of  the  initial  relative  humidity  of  the  clearing,  for 
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various  values  of  C - C,  and  T - T,.  The  lower  diagram  shows  the  same  for 
O 1 O 1 

a slot,  or  swath,  of  clearing,  which  was  the  particular  geometric  shape  discussed 

1 5 2 

by  Elliott.  The  K value  assumed  in  the  plotting  was  10  cm  /sec.  For  K 
e ^ 

values  other  than  this,  it  might  be  mentioned  that  the  closing  times  vary  inversely 

with  K . 
e 

The  diagrams  reveal  that  the  closing  times  are  sensitively  dependent  on  the 

I initial  relative  humidity  of  the  clearing,  particularly  in  the  humidity  range  from 

90  to  100  percent.  They  also  demonstrate  that  a circular  clearing  — for  a circle 
diameter  equal  to  the  slot  width  — will  close  about  twice  as  rapidly  as  a slot  of 
clearing. 
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CYLINDER 


Relative  Humidity— »60% 


t.  minutes 


EXISTENCE 
OF  CLEARING 
IMPOSSIBLE 


SLOT 


Initial 

Relative  Humidity— •60% 


t.  minutes 


EXISTENCE 
OF  CLEARING 
IMPOSSIBLE 


Initial  Woter  - Substance  Difference,  C0-C,  , gm/m3 


Figure  3.  Closing  Times  for  a Cylindrical  Clearing  (upper 
diagram)  and  for  a Slot,  or  Swath,  of  Clearing  (lower  dia- 
gram). The  diagrams  show  how  the  closing  times  vary  as 
a function  of  the  initial  relative  humidity  for  particular 
values  of  CQ  - Cj  and  T0  - Tj.  The  diagrams  are  drawn  for 
Ke  = 10'  cm^/sec,  Tj  ; 10°C,  and  Cj  9.6  gm /m^,  with 
Mj  = 0.2  gm/m  . The  initial  diameter  of  the  cylinder  was 
assumed  to  be  200  meters,  which  also  corresponds  to  the 
assumed  initial  width  of  the  slot.  The  Sc  values  indicated 
are  values  or  (C  - C )/(CQ  - C^  (TQ  - TJ/(TQ  - Tt). 

The  A,  B,  and  C points  shown  in  the  upper  diagram  corre- 
spond to  the  "warm  dry,"  "warm  moist,"  and  "cold  dry" 
situations  illustrated  in  Figure  2. 
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